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ABSTRACT 

The  general  dif ferential  equations  describing  the  behavior  of  the 
rotary  regenerator  including  longitudinal  heat  conduction  in  the  direction 
of  fluid  flow  are  sufficiently  complicated  to  preclude  a  complete  analyti- 
cal solution.  A  few  solutions  found  in  the  literature  restricted  to  very 
special  cases  are  discussed. 

A  numerical  finite  difference  method  is  presented  which  will 
determine  the  effect  of  longitudinal  heat  conduction  in  rotary  regenerators 
for  steady  state  conditions.   In  the  development  no  assumptions  are  made 
which  would  restrict  the  range  of  parameters  for  which  the  analysis 
would  be  applicable. 

The  conduction  effect  on  the  regenerator  effectiveness  was  evaluated 
with  the  view  of  obtaining  results  most  useful  for  the  gas  turbine 
regenerator  problem,  however,  these  results  may  also  be  used  for  other 
regenerator  problems. 

A  CDC  1604  digital  computer  was  used  to  carry  out  the  computations. 
The  results  are  presented  graphically  and  in  tabular  form,  employing  a 
suitable  set  of  non  dimensional  parameters.   The  range  of  parameters 
which  have  been  covered  are: 

0.5  4  Cmin/Cmax^l.O 
1-0   Gr/CmtnaO.O 
1.0   NTUo-20.0 
0.25  4  (hA)*  1  ..: 
0        0.2 
0.15  As*  1.3 
The  author  expresses  his  appreciation  to  C.  P.  Howard,  Associate 
Professor,  for  his  direction  and  encouragement  in  this  work. 
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of  Qaax.  »q  ft 
Asn  -  solid  area  available  for  longitudinal  heat  conduction  on  the  side 

of  C,nin,  sq  ft 
As  =  total  solid  area  available  for  longitudinal  heat  conduction,  sq  ft 
C  *  heat  capacity  rate  (Wc)  of  fluid  or  rotor  matrix  according  to 
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c  ■  specific  heat  of  fluid  (constant  pressure)  or  rotor  matrix  material 

depending  on  subscript,  Btu/lb,  deg  F 
h  ■  unit  conductance  for  thermal  convection  heat  transfer,  Btu/hr, 

deg  F,  sq  ft 
k  ■  unit  thermal  conductivity,  Btu/hr,  aq   ft,  deg  F/ft 
L  ■  total  length  of  the  matrix  in  the  direction  of  fluid  flow,  ft 
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T  =  temperature  of  matrix,  deg  F 

l'x(i,j)  =  subscripted  matrix  temperature  on  the  sldr.  of  C^*,  deg  F 
T  . f  .  ■  subscripted  matrix  temperature  on  the  side  of  C^in*  deg  F 
t  ■  temperature  of  fluid,  deg  F 

tx(i  j\  »  subscripted  fluid  temperature  on  the  side  of  OWx,  deg  F 
*n(f  a)  "  subscripted  fluid  temperature  on  the  side  of  C^n,  deg  F 
W  =  mass  flow  rate  of  fluid  (lb/hr)  or  matrix  (lb,  rev/hr),  according 

to  subscript 
y  ■  distance  of  a  point  from  the  cold  end,  measured  in  units  of  total 

length  of  matrix,  0  y  L 
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I.    Introduction 

The  regenerative  type  of  heat  exchanger  has  been  subjected  to 
extensive  treatment  within  the  last  thirty  years  or  so  and  contributions 
to  the  theory  of  regenerators  have  been  made  by  many  writers.   The  most 
questionable  of  the  idealizations  made  in  the  development  of  the  theory 
is  that  of  zero  matrix  conductivity  in  the  direction  of  fluid  flow.   The 
loss  due  to  heat  conduction  in  the  matrix  material  is  specially  important 
in  the  case  of  high  performance  regenerators  in  which  a  small  reduction 
in  the  effectiveness  will  result  in  a  large  increase  of  the  heat  transfer 
areas  to  compensate  for  the  loss.  Another  important  factor  is  the  modern 
trend  toward  smaller  sizes,  since  the  loss  due  to  conduction  of  heat  in 
the  matrix  obviously  increases  with  decreasing  length. 

Schultz   2  '  solved  the  problem  including  conduction  for  the 
special  case  of  the  balanced  regenerator  and  very  large  values  of  Cr/C^n 
(i.e.,  the  capacity  rate  ratio  of  rotor  matrix  to  minimum  fluid  capacity 
rate).   In  this  case  the  system  of  partial  differential  equations  defining 
the  problem  for  the  steady  state  are  reduced  to  ordinary  differential 
equations  and  the  problem  can  be  solved  analytically. 

Hahnemann  3  extended  the  solution  to  the  unbalanced  regenerator 
(i.e.,  unequal  partitions  of  the  cross  sections  and  unequal  heat  transfer 
coefficients  together  with  unequal  capacity  rate  ratios  of  the  fluid 
streams)  but  still  for  an  infinitely  large  value  of  Cr/C^^. 

London  has  shown,  as  reported  by  Lambertson  1   ,  that  a  correction 
factor  for  conductivity  in  the  direction  of  fluid  flow  for  approximately 

^Number  in  brackets  refers  to  bibliography  on  page  2.1 


equal    fluid  capacities,    of  the    form: 

Ae  _   I 


C  L 
should  result  in  a  somewhat  pessimistic  prediction  of  the  reduction  in 
effectiveness. 

It  is  the  objective  of  this  thesis  to  develop  and  present  a  finite 
difference  numerical  analysis  with  solutions  to  solve  for  the  conduction 
effect  on  effectiveness  for  the  rotary  regenerator.   Since  there  are  six 
dimensionless  parameters  required  to  specify  behavior  in  the  most  general 
case,  the  ranges  of  the  parameters  covered  in  the  solutions  will  be 
those  of  most  interest  to  the  gas  turbine  regenerator,  however,  the 
analysis  will  place  no  restrictions  on  these  parameters. 


2.   Method 

The  conventional  idealizations  and  boundary  conditions  assumed  in 
the  derivation  of  the  governing  differential  equations,  including  the 
heat  conduction  in  the  direction  of  fluid  flow  4   ,  are  the  following: 
(See  Fig.  la) 

a.  The  thermal  conductivity  of  the  matrix  material  is  zero  in  the 
direction  of  matrix  metal  flow.   It  has  a  finite  value  in  the 
direction  of  fluid  flow,  and  is  infinite  in  the  other  direction 
normal  to  the  fluid  flow. 

b.  The  specific  heats  of  the  two  fluids  and  matrix  material  are 
constant  with  temperature. 

c.  No  leakage  of  the  fluids  occurs  either  due  to  direct  leakage  or 
carry  over,  and  each  fluid  flow  is  unmixed. 

d.  The  convective  conductance  between  the  fluid  and  the  matrix  are 
constant  with  flow  length. 

e.  The  fluids  pass  in  counter  flow  directions. 

f.  Entering  fluid  temperatures  cire  uniform  over  the  flow  inlet  cross 
section  and  constant  with  time. 

g.  The  matrix  temperature  gradient  in  the  direction  of  fluid  flow  is 
zero  at  the  ends. 

h.  Regular  periodic  conditions  are  established  for  all  matrix 
elements,  i.e.,  steady  state  condition. 

The  differential  equations  derived  by  Schultz  2  are  given  in 
Appendix  1 . 

It  was  suggested  b>  Do  ..'   rre  *nd  carried  out  by  Lambert  son  [1  i  for 
the  case  of  no  conduction,  that  a  rotary  regenerator  could  be  divided 


Typical 
element 


FIG.  la   ILLUSTRATIVE  MATRIX  ARRANGEMENT 
AND  FLUID  FLOW 


Lnto  finite  elements  a^  represented  schematically  In  Fig.  lb.  This  same 
method  c/in  be  extended  to  the  case  with  conduction,   if  t'.io  elements  are 
considered  to  be  fixed  in  space,  each  element  can  be  regarded  as  a  cross 
flow  heat  exchanger  with  a  ges   ytream  and  a  metal  stream.   For  an  element 
on  the  side  of  C^^,  the  heat  transfer  rate  by  convection  is  equal  to  the 
rate  of  change  in  enthalpy  of  the  fluid  across  the  element, 


AT     -    I     -T 


(1) 


where  _T    represents  the  mean  temperature  difference  between  the  fluid 
avg 

and  the  matrix  element.  For  a  siall  enough  element  the  arithmetic  mean 
temperature  difference  may  be  assumed  valid,  so  that 


A  T     =  —  +  t 


X(x 


where  the  fluid  and  matrix  tempe-atures  are  subscripted  to  indicate  the 

average  temperatures  across  the  inlet  and  outlet  of  the  element  for 
simple  unmixed  cross  flow. 

Considering  an  energy  balance  for  the  element,  the  energy  transferred 

to  the  element  by  convection  plus  energy  transferred  to  the  element  by 
conduction  less  the  energy  transferred  out  of  the  element  by  conduction 
must  equal  the  energy  stored  in  the  element: 


^, 


. 


tx     -t 


H-  I 


all :. 


T 

i1  +  ■ 


r 

.  ,  '  >  L-  ■'  :^'   ^'  <"'j  ~~  I   J  _  Cr 


,1 


Mr 


Nr 


(2> 


By  using  equations  (1)  and  (2)  the  outlet  temperatures  for  an 

clement  on  the  side  of  C^y   laay  be  solved  for  in  terms  of  the  remaining 
temperatures  and  the  result  expressed  in  dimensionless  form  (Appendix  1) . 

Since  the  longitudinal  heat  conduction  is  zero  at  the  ends,  the 
energy  balance  equation  is  different  from  equation  (2)  for  the  elements 
of  the  first  and  last  row.   Three  different  expressions  for  the  matrix 
outlet  temperatures  are  obtained: 
First  row  elements 


|)    I     !)  +R?  I  «)  J 


(3) 


Middle  row  elements 
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Lsst  row  elements 
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By  solving  for  the  fluid  outlet  temperature  on  the  si.de  of  Cmax  in  terms 
of  the  matrix  outlet  temperature  from  equation  (1),  only  one  expression 


is  obtained,  which  holds  for  all  the  elements  on  the  side  of  C 


max 


*■  l> 
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Similarly,    tor  c*  element  on  the  side  of  C,nin 
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with 
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Energy  balance. 
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By  the  same  method  the  outlet  temperatures  of  the  elements  on  the  side 
of  C^j,  may  be  obtained  using  equations  (7)  and  (P.). 
First  row  elements 


Middle  row  elements 


(9) 


=  F  f    -  F  \~ 
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(10) 


L&st  row  elements 


' 


(H) 


Solving  for  the  fluid  outlet  temperature  for  the  elements  on  the  side 
of  Cfnin   from  equation  (7)  gives 


f    3   ~  6t^.;'  +  ^Umm)  +  ^.M+i 


(12) 


■CF/3) 
which  holds  for  all  the  elements  on  the  side  of  C^j n.   For  the  values 

of  the  constants  see  Appendix  1. 

In  addition  to  the  four  dimensionless  groupings  for  no  conduction, 

Gmin/Cmax  =  capacity  rate  ratio  of  fluid  flow  streams 

Cr/Cm^n  =  capacity  rate  ratio  of  rotor  matrix  to  minimum 

fluid  capacity  rate 

* 

(hA)   ■  conductance  ratio 

NTUo  ■  ,  over-all  number  of  transfer  units; 

two  additional  dimensionless  groupings  are  needed  in  the  coefficients 
of  equations  (3),  (4),  (5),  (6),  (9),  (10),  (11),  and  (12)  to  take 
conduction  into  account.   These  conduction  parameters  may  be  defined  as: 


,  conduction  parameter  on  side  of  Cmax 

i_ 

. 

,  conduction  parameter  on  side  of  C^,in. 


Also  


since  the  thermal  conductivity  and  the  length  of  the  matrix  are  the  sa^-.e 
for  both  sides,  and  where 


,  conduction  area  ratio. 
Another  form  of  defining  the  conduction  parameter  is 


-  -  ,  total  conduction  parameter 

In  this  case  As  is  the  total  solid  area  available  for  conduction,  that  is 

This  form  offers  the  advantage  that  for  a  given  value  of  the  conduction 
parameter,  the  resulting  effectiveness  will  not  be  affected  by  small 
changes  in  As  . 

From  the  schematic  representation  shown  in  Fig.  lb,  it  should  be 
noted  that  the  left  edge  of  the  regenerator  is  physically  the  same  as 
the  right  edge  and,  therefore,  the  matrix  inlet  temperature  for  the 
elements  of  the  first  column  on  the  side  of  C^x  is  th^  same  as  the 
outlet  temperature  of  the  corresponding  element  of  the  last  column  on 
the  side  of  C min*   Expressed  mathematically, 

= 

where  1  and  f  can  take  the  values  of  1,  2,  3, Kr, 

This  is  referred  to  os  the  reversal  condition. 

From  equations  (3),  (4),  (5),  (9),  (10),  and  (11)  it  is  seen  that 
in  order  to  solve  for  the  outlet  matrix  temperature  for  a  particular 
element,  it  is  necessary  to  know  the  matrix  temperatures  of  the  next 
element.  An  estimate  of  those  values  can  be  made  by  determining  the 


temperature  distribution  for  the  cas*:  of  aero  conduction  In  the  (Ui.fr 
of  fluid  Clow,  since  the  constants  D3  and  Dg  In  equations  (3),  (4)  and 
(5)  and  the  constants  F3  and  Fg  In  equations  (9),  (10)  and  (11)  are  zero 
in  this  case. 

A  temperature  scale  can  be  used  for  which  the  fluid  entrance  tempera- 
ture is  zero  at  Che  side  of  CtlIn  and  unity  at  the  side  of  C^x' 

tn(l,g)  =  °     for  8-1.2,3, Nn 

^(ij)  "  l      for  j  =  1,2,3, Nx 


For  calculation  purposes  a  temperature  distribution  is  assumed  on 
the  left  edge  and  the  problem  is  solved  for  the  no  conduction  case  in 
order  to  obtain  an  initial  estimate  of  matrix  temperatures  for  working 
the  problem  with  conduction.  This  also  provides  the  no  conduction 
effectiveness  necessary  for  comparison.. 

The  outlet  temperatures  are  calculated  for  every  element  by  repetitive 
use  of  equations  (3),  (A),  (5)  and  (6)  depending  on  the  location  of  the 
element.  The  calculation  is  started  with  the  first  element,  first 
column  on  the  side  of  C^ax  and  then  xx>rking  down  the  column.  When  the 
first  column  of  elements  is  completed,  the  second  column  is  calculated 
and  so  on  until  the  outlet  temperatures  of  all  the  elements  on  the  side 
of  Cmax  are  determined.   It  was  found  that  convergence  could  be  enhanced 
if  two  passes  per  column  were  made  before  proceeding  to  the  next  one. 
The  reason  for  this  is  that  the  equations  for  the  outlet  temperatures 
contain  the  matrix  temperatures  of  the  next  element  which  are  necessary 
to  estimate. 

At  the  seal  represented  by  the  double  line  in  Fig.  lb, 

I 


Tx(i,Nx  +  1)  "  Tn(f,l) 
since  the  matrix  outlet  temperature  for  the  elements  of  the  last  column 
on  the  side  of  Qnax  *-8  physically  the  same  as  the  inlet  temperature  of 
the  corresponding  element  of  the  first  column  on  the  side  of  C^^. 

The  same  method  is  applied  to  the  side  of  Cg^n  as  to  the  side  of 
CMax.   If  the  temperature  distribution  assumed  on  the  left  edge  was 
correct  it  would  then  be  duplicated  on  the  right  (i.e.,  the  reversal 
condition  is  fulfilled).   If,  however,  this  is  not  the  case,  the 
resulting  temperatures  on  the  right  side  Tn(f^«  +j\ |  are  now  used  on 
the  left  ,Txm  j\J  and  the  procedure  is  repeated.  After  each  pass  (i.e., 
the  complete  calculation  of  a  temperature  distribution)  an  energy 
balance  is  made;  and  before  the  solution  is  accepted  for  a  particular 
set  of  parameters,  the  heat  balance  error  together  with  the  reversal 
condition  has  to  be  fulfilled  to  the  specified  accuracy. 

The  effectiveness  of  a  heat  exchanger  is  defined  (Appendix  2)  as 


= 


C-nwM.  ^Txi  —  X/rui.  I  ^X*    1*1 

From  the  conditions  of  the  problem  the  above  expression  for  effectiveneso 
reduces  to 


e  - 


r 

.Nr+l, 


3) 


a1 
The  heat  balance  error  was  computed  from 


. 


!  -[Wwi) 
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That  1    solution  will        -  '  I  .   ,  .  »*■  r-  •*•?  1  ,  i  I 

•oi.t  »u lance  error  converge  to  zero)  will  depend  on  t'rns  nuuiber 
of  elements  use-?.  More  important  than  this  is  the  relation  existing 
between  the  number  of  subdivisions  for  the  three  streams.  The  sufficient 
but  not  nec<?ssnr>  conditions  for  convergence  are  (Appendix  1)  that 

/       i 

M   >(A^  J 

\As*J 

/^  ■ 

Fig.  2  is  a  plot  of  the  different  values  of  effectiveness  obtained 
for  the  following  set  of  parameters- 

Several  curves  are  obtained  depending  on  the  relation  existing  between 
the  number  of  subdivisions  for  every  stream.   To  extrapolate  the 
effectivewass  to  an  infinite  number  of  elements,  these  curves  no.  » e 
approximated  by  straight  lint's.   The  error  introduced  will  depend  on 
the  number  of  elements  usod  and  will  not  be  constant  for  all  th»?  range 
of  the  partneters.   As  on  illustration  consider  the  set  of  parameters 
used  to  plot  Fig.  2.  Three  cases  are  investigated: 

1)  =  Nn  -  «r/2 

2)  -  Nn  •  Nr 

3)  L  «  Nn  -  2Hr 

A  linear  extrapolation  is  performed  using  the  values  of  effectiveness 

obtained  for  values  of .   From  the  results  it  is  founJ 

Nr(Nx  +  Nn) 

that  the  minimum  Mr  necessary  to  give  an  error  in  the  extropoiated 
<  tectiveness  of  less  than  0.01%  is  18,  for  which   ev .  -  76.57X 


1% 


%  SS3N3AU.033J3 


u 


lor  the  three  cas  s.   L'uc  cioLoi'.flinlnj*  factor  in  the  extrapolation  is 
the  nunber  of  matrix  elements.   The  higher  the  i/aluc  oi  Nr  used  for  the 
extrapolation  the  smaller  the  error  in  the  straight  line  approxii.<atio n 
Unfortunately,  an  increase  in  Nr  will  bring  a  higher  increase  of  1>X   anu 
Nn,  since  the  convergence  condition  for  this  example  is 

Nx  =  Mn  >  o.of  N/  . 

That  this  is  only  a  sufficient  condition  can  be  seen  from  the  case  of 

Nx  -  Nn  -  Nr/2 
for  which  the  solution  converges  for  values  of  Nr  5;  and  when 

Nx  -  Nn  -  Nr 

the  solution  will  still  converge  for  Hr  20.   It  is  true,  however,  that 
the  closer  the  number  of  matrix  elements  are  to  the  sufficient  convergence 
condition,  the  fewer  passes  it  takes  for  the  solution  to  converge  to 
the  specified  error  in  the  effectiveness.   Table  1  shovs  the  number 
of  passes  necessary  to  obtain  a  ,  olution  with  a  heat  balance  error  of 
less  than  0.0057.. 

For  a  given  set  of  parameters  the  effectiveness  was  determined 
first  using  one  set  of  subdivisions  and  then  recalculated  increasing 
the  number  of  subdivisions.  A  linear  exirapolation  was  performed  to  obtain 
the  effectiveness  corresponding  to  an  infinite  number  of  elements.   The 
accuracy  for  each  value  of  effectiveness  was  specified  to  be  of  four 
significant  figures.   When  these  values  were  extrapolated  to  an  infinite 
number  of  elements,  the  accuracy  was  reduced  to  only  three-  significant 
figures  since  a  compromise  had  to  be  riade  with  the  time  available  for 
computations  and  the  number  of  runs  to  be  made.   Howev?r,  ti  *  iet?iod 

L5 
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Number  of  passes  necessary  to  obtain  an  effectiveness  with  a  heat 
balance  error  of  less  than  0.005%  for  the  follo;:i.n°  stt  of  parameters: 
Cbln/Cmax  "1*0  ,  ^/C^,  =  2.0  ,  <hA)*  =  1.0  , 
trruo  =  10.0  ,     \«  0.1  ,   As*  =  1.0 


Nr 

Nx 

Nn 

-  ■  ' ■  -   '       .........  . —    ._| 

Mo.  of  passes 

* 

e  .   /» 

12 

12 

12 

27 

76.581 

12 

24 

24 

22 

76.581 

14 

7 

7 

41 

76.578 

14 

14 

U 

34 

76.576 

14 

28 

28 

27 

76.576 

16 

16 

16 

41 

76.574 

16 

32 

32 

32 

76.574 

18 

9 

9 

62 

76.573 

18 

18 

13 

50 

76.572 

ia 

36 

3C 

38 

76.572 

20 

20 

20 

60 

76.572 

22 

22 

22 

71 

76.571 

*  e^  is  obtained  by  a  linear  extrapolation  of  the  effectiveness 


for  Nr  and  Mr-1 
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oses  o  cestri  stioti  as  to  Che  ul     tc  accuracy  of  the  solution. 

When  working  with  very  high  values  of  Cr/C,n£n  the  computation  time 
becomes  considerably  long.   This  is  due  to  the  fact  that  the  final 
temperatures  differ  very  little  from  the  assumed  initial  temperatures  and 
a  large  number  of  elements  is  required.  The  starting  estimates  of  matrix 
temperatures  then  become  very  important.   For  the  limiting  case  of  Cy/C^n 
equal  to  infinite,  there  is  no  temperature  variation  in  the  direction 
of  metal  flow  so  that  the  method  without  modifications  is  not  applicable. 
For  this  case,  however,  the  exact  solution  of  Hahnemann  [3 J  is  available. 

The  Fortran  system  i«is  used  to  program  the  above  method  for  the 
CDC  1604  digital  computer.  The  matrix  notation  employed  in  the  equations 
derived  previously  is  very  convenient  in  writing  up  the  program,  since 
several  loops  can  be  made  by  just  changing  the  subscripts  of  the  variables 
involved.  The  complete  computer  program  is  given  in  Appendix  2. 
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3.      us  8  ion    'esults 

The  solutions  generated  are  presented  tn  tabular  fon:i  Ln  Tables  2 
through  25  and  graphically  in  Figs.  3  and  4.  The  accuracy  in  the 
effectiveness  was  specified  to  be  of  three  significant  figures.   It  was 
found  that  for 

0.7  (hA)*\  As*  1.3  (hA)* 
the  resulting  values  of  effectiveness  were  almost  identical  to  the  ones 
calculated  for 

As*  -   (hA)* 
so  that  these  were  the  values  used  in  preparing  the  Tables. 

By  using  NTUo  and  A  as  defined  in  the  nomenclature,  which  contain 
(hA)  and  As'  explicitly,  the  influence  of  these  two  parameters  on  the 
conduction  effect,  -3-  ,  were  found  to  be  very  small  for  values  of 
^min/^taax   ^.9  and  \  0.1,  which  is  the  range  for  the  gas  turbine 
regenerator  problem. 

Figs.  3  and  4  show  the  effect  of  conduction  on  effectiveness  for 
values  of  Cn^n/C^ax  =  1.0  and  0.9  respectively,  with  different  values  of 
NTUo,  Cf/Cj^n  and   ,  up  to  a  limit  of  107..   The  conduction  effect 
continues  to  increase  with  increasing  A •   For  small  values  of  A  this 
effect  increases  nearly  exponentially.  By  further  increasing  A  .  a 
saturation  is  reached  which  approaches  the  limiting  condition  of  \  =  oo. 
Fig.  5  is  a  representative  curve  which  shows  the  limiting  influence  on 
the  conduction  effect  by  continuously  increasing  \  for  a  given  set  of 
parameters. 

The  conduction  effect  for  values  of  Cr/Cm^n  =  10. 0  are  almost  identical 
with  the  values  reported  by  Schultz  [2  J  and  Hahnemann   [3  J  ,  [6  J  for  the 
limiting  case  of  Cy/C^n  =  oo. 
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For  small  values  of  NTUo  the  effect  of  increasing  Cr/Cain  la  to 
increase  slightly  the  conduction  effect,  as  can  be  seen  from  Fig.  6. 
For  higher  values  of  NTUo  the  effect  of  conduction  is  reduced  with 
increasing  ^/C^^. 

The  effect  of  increasing  NTUo  upon  ^ is  intimately  related  to  the 
value  of  Cr/Cn,in  and  Cmln/Cmax.   This  effect  can  best  be  explained  by 
observing  Fig.  6.   For  the  case  where  Cxnin/Cjaa^   =  1.0  the  effect  of 
increasing  NTUo  will  result  in  an  increase  in  4£  regardless  of  the  value 
of  Cr/C,nln;  of  course,  the  higher  the  value  of  Cr/Cmin  the  smaller  this 
effect  will  be  for  relatively  high  values  of  NTUo.   Now,  by  reducing 
Cm£n/Cmax  and  increasing  Cr  /Gn,£n  a  point  is  reached  in  which  an  increase 
in  the  NTUo  will  result  in  a  reduction  in  ^  .   The  point  where  this 
occurs  will  depend  upon  the  value  of  the  conduction  parameter  X  . 
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4.   Summai         •  is  ions 

A  finite  dis     ce  numerical  analysis  Is  presented  which  will 
*;.  termine  the  effect  of  the  longitudinal  heat  conduction  in  rotary 
regenerators  for  steady  state  conditions. 

The  ranges  of  governing  parameters  which  have  been  covered  are: 

0-5  ^W(W  41.0 

l-°     cr/cmin^°-° 
1.0  iUTUo  <20.0 

0.25  4 (hA)*4  1.0 

0  ^  A  ^0.2 

O.lSiAs*  ^1,3 

The  ranges  of  parameters  of  interest  for  gas  turbine  regenerator 
application  are  plotted  in  Figs.  3  and  4.  A  more  accurate  determination 
of  the  conduction  effect  than  that  given  in  the  figures  can  be  made  by 
interpolating  values  from  the  Tables. 

The  computer  program  used  to  carry  out  the  calculations  on  the 
CUC  1604  digital  computer  is  included  in  Appendix  2.   This  program 
imposes  no  restrictions  on  the  p..tasiieters  for  which  the  analysis  would  be 
applicable.   It  is  recommended,  however,  that  for  values  of  Cr/Cmin>10 
the  solution  for  Cj. /Cmin   "  10  be  used.   For  this  case  the  computer  time 
for  a  solution  increases  considerably  and  the  error  introduced  by  assuming 
Cr/Cmin  =  l0  is  negligible. 
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TABLE   2  I'UDTNAL  MEAT  CONDUCTION  EFFECT 


C         / 1 ' 
"V.in     'max 


1.0    ,    Cj/CiajLn  -  1.0    ,    (hA)      -  As*  -  1.0 


e  % 


Conduction  effect,---    % 


NTUo 


,A 


0.01       0.02 


0.04 


0.08 


0.14 


0.16 


0.20 


4 
5 
6 
7 
8 
9 
10 

12 
14 
16 
18 
20 
40 


46.64    i  0.42 


60.06 


3     :   66.72 


0.58 
0.68 


70.86    ;  0.77 

73.75    |  0.84 

i 

I 

75.92 
77.63 
79.01 
80.16 

81.15    |  1.12 

i 
82.74 

84.00 

85.01 

85.86 

86.58    i   1.49 

90.37       1.92 


0.78 
1.11 
1.32 
1.48 
1.63 
1.76 
1.87 
1.97 
2.07 
2.16 

2.32 
2.46 
2.60 
2.72 
2.83 
3.55 


1.48 

2.50 
2.82 
3.09 


4.05 


5.17 
6.30 


2.61 
3.81 
4.58 
5.17 
5.65 
6.05 
6.41 
6.72 
7.01 
7.26 

7.71 
8.10 
8.44 
8.74 
9.01 
10.54 


3.94 

5.93 

7.17 

8.08 

8.81 

9.42 

9.93 

10.38 

10.78 

11.13 

11.74 
12.25 
12.69 
13.08 
13.42 


4.30 
6.54 
7.92 
8.93 
9.73 


12.23 


14.60 
16.58 


4.95 
7.65 
9.30 
10.48 
11.41 
12.16 
12.80 
13.34 
13.82 
14.24 

14.95 
15.54 
16.04 
16.47 
16.85 


0.32 


6.42 
10.03 
12.60 
14.22 
15.44 


18.97 


21.96 
24.25 


100 

!  93.87 

2.70 

4.78 

7.97 

12.55 

18.79 

26.54 

500 

97.87 

4.04 

6.80 

10.57 

15.42 

21.65 

29.27 

1000 

98.81 

4.39 

7.30 

11.17 

16.07 

22.30 

29.88 

TABLE   3 

LONGITUDINAL  HEAT  CONDUCTION  EFFECT 

1 

Tniir  *tnax 

■U,   CrAWn  «  1-2   ,    (hA)*  «  A»*  «  1 

.0 

e 

Conduction  effect    —  % 

NTUcX  A 

0 

0.02                 0.08                 0.14 

0.20 

1 

47.62 

0.81                  2.66                 4.01 

5.04 

2 

61.97 

1.13                  3.88                  6.05 

7.81 

3 

69.15 

1.33                 4.64                  7.28 

9.46 

4 

73.61 

1.49                  5.21                  8.18 

10.63 

5 

76.72 

1.62                  5.66                  8.89 

11 .  54 

6 

79.04 

1.73                  6.05                  9.47 

12.28 

7 

80.85 

1.34                  6.38                  9.96 

12.89 

8 

82.31 

1.93                 6.67               10.38 

13.40 

9 

83.53 

2.01                  6.03                10.76 

13.86 

10 

84.55 

2.09                  7.1c                11.08 

14.26 

11 

85.44 

2.16                  7.37                11.38 

14.61 

12 

86.21 

2.22                 7.56               11.65 

14.93 

13 

86.89 

2.28                  7.74                11.89 

15.22 

14 

87.49 

2.34                  7.90                12.11 

15.48 

15 

88.03 

2.40                  8.05                12.31 

15.71 

16 

88.52 

2.45                  8.19                12.50 

15.93 

17 

83.96 

2.50                  3.32                12.67 

16.13 

18 

89.37 

2.54                 8.44               12.83 

16.31 

19 

89 .  74 

2.59                  8.56               13.00 

1 6 .  48 

20 

90.09 

2.63                  3.66                13.12 

16.64 

4     LONGITUDINAL  HEAT  CONDUCTION   EFFECT 
<Wa'<W  "   l'°    '    Cr/Cmin  "   1 '4    ■    0*>*  =  As*  =   1 


e  % 

Conduction  effect  ( 

NTUo' 

1 

0 

.01    .02   .04    .08 

.14  .1 

16    .2    .32 

|    1 

48.23 

.82       2.69 

4.05 

5.10 

2 

63.16 

1.14       3.93 

6.12 

7.90 

3 

70,65 

1.33       4.67 

7.33 

9.  S3 

!     * 

75.27 

1.48       5.20 

8.20 

10.67 

5 

78.45 

1.60       5.63 

8.86 

11.54 

6 

80.82 

1.70       5.98 

9.40 

12.23 

7 

82.65 

1.79       6.28 

9.85 

12.80 

8 

84.12 

1.87       6.54 

10.24 

13.28 

9 

85.33 

1.94       6.76 

10.57 

13.69 

10 

86.34 

2.00       6.96 

10.86 

14.05 

12 

87. 9G 

2.12       7.31 

11.36 

14.64 

14 

89.20 

2.22       7.59 

11.75 

15.12 

16 

90.18 

2.3        7.83 

12.08 

15.51 

18 

90.98 

2.37       8.04 

12.36 

15.84 

20 

91.65 

2.44       8.22 

12.61 

16.12 

24 


TABLE  5  LONGITUDINAL  HEAT  CONDUCTION  EFFECT 


Cmln/Cmax  -  1.0   ,   Cr/C^m  -  2.0   ,    (hA)*  =  As*  •  1.0 


e  % 


Conduction  effect,  ~  % 


NTVo 


,01 


.02 


,04 


.08 


.14 


16 


.2 


.32 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

12 

14 

16 

18 

20 


49.12 
64.91 
72.80 
77.60 
80.86 
83.22 
85.03 
86.46 
87.63 
88.59 
90.11 
91.25 
92.14 
92.86 
93.45 


500         99 . 65 


lono 


99.80 


.84 
1.16 
1.34 
1.47 
1.56 
1,64 
1.71 
1.76 
1.82 
1.86 
1.94 
2.01 
2.06 
2.11 
2.15 

1.40       2.72 

1.40       2.72 


5.11 


5.13 


2.73 

4.11 

4.0 

6.22 

4.70 

7.40 

5.18 

8.19 

5.54 

8.78 

5.83 

9.23 

6.07 

9.61 

6.27 

9.92 

6.44 

10.18 

6.60 

10.41 

6.85 

10.79 

7.05 

11.08 

7.22 

11.33 

7.36 

11.53 

7.49 

11.71 

9.10 


9.14 


15.18 


15.24 


5.17 
8.04 
9.37 
10.70 
11.47 
12.07 
12.55 
12.95 
13.28 
13.57 
14.04 
14.41 
14.71 
14.95 
15.16 


23.40 


23.47 
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TABLE   6     "  IDINAL    IEAT  CONDUCTION   ' FFECT 

Cmin/Cmax  -  1.0  ,   Cr.'Cmin  -  10.0   ,    (hA)*  -  As*  =  1.0 


e  % 


Conduction  effect  -  °7 


NTUo 


.02 


.08 


14 


.2 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 

12 

14 
16 
18 
20 


49.96 
66.59 
74.91 
70.89 
83.22 
85.59 
87.38 
88.76 
89.87 
90.78 

92.18 
93.20 
93.99 
94.61 
95.11 


.85 
1.18 
1.35 
1.46 
1.53 
!  .58 
1.62 
1.65 
1.68 
1.70 

1.73 

1.76 
1.76 
1.80 
1.81 


2.77 
4.04 
4.72 
5.13 
5.42 
5.62 
5.78 
5.91 
6.02 
6.10 

6.24 
6.34 
6.43 
6.49 
6.55 


4.16 
6.29 
7.41 
8.10 
8.57 
8.92 
9.19 
9.40 

57 

9.71 

9.94 
10.11 
10.25 
10.36 
10.45 


5.24 
8.12 
9.63 
10.57 
11.21 
11.67 
12.03 
12.31 
12.54 
12.73 

13.03 
13.26 
13.45 
13.59 
13.71 
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BLE  7     LONGITUD]    " 

:;.-.>•     IONDUCT 

e  % 

tax  "          ' 

cr /Qnin 

«  1.0   ,    (hA)*  = 

«o.5 

Conduction  effe  t        -  ^zl 

_ 

NTUo 

0 

0.015 

.03 

.06          o.:^ 

0. 

0.*8 

■                                     1 

1 

46.62 

.62 

1.17 

2.12             3.64 

5.72 

8.05 

2 

60.08 

.85 

1.64 

3.04              5.36 

.03 

3 

66.75 

1. 01 

1.95 

3.63             6.46 

10.67 

]  5 .  90 

4 

70.90 

I .  L4 

2.19 

4.09              7.26 

11. 

1   1      p  '■ 

L  7 .  < 

5 

73.79 

1.25 

2.40 

4.47              7.91 

12.98 

19.26 

10 

81.19 

1.66 

3.15 

5.77              9.98 

15. 

21 

20 

86.62 

2.19 

4.07 

7.22            12. 

IS. 

26, 

40 

90.46 

2.69 

4.92 

8.51           13.77 

l^.Tj 

• 

28  J 

100 

93.61 

3.24 

5.34 

9.85           15.4 

22.60 

64 

500 

95.57 

3.51 

6.36 

10.66 

23.  "r 

31. .  91 

iooo 

95.83 

3.55 

6.43 

10.77 

92 

32.10 

27 


• 

Cmin/Cmax  -  1.0  ,   Cx/Cmiri  :    2 

0    , 

i 

.5 

!      e  % 

i              . 

Conduction 

ect 

> 

KTUo\ 

\ 

0 

.015          .03            .06 

.12 

.21 

1 

49,11 

1.22 

3.79 

5.51 

6.75 

2 

64.92 

1.70 

5,62 

8.47 

10 .  64 

3 

72.81 

1.97 

0.:-', 

10.10 

12.78 

4 

77.61 

2.16 

7.32 

.17 

\L.\1 

5 

80.86 

2.30 

7.02 

ll  .9: 

15.1.6 

6 

83.23 

2.42 

8.22 

12.55 

L5.92 

7 

85.04 

2.51 

8.54 

13.03 

16.52 

8 

86.47 

2.60 

8.81 

13. 

17.01 

9 

87.63 

2.67 

9.04 

13 

1 7  .  '.  2 

10 

88.59 

2.73 

9 .  24 

14.05 

I  ) 

12 

90.11 

2.84 

9.57 

14.52 

L8.32 

14 

91.25 

2.94 

9.83 

14.88 

L8.76 

16 

92.14 

3.01 

10.05 

15J 

18 

92.86 

3.08 

10.23 

15. 

19.40 

20 

93.45 

3.14 

10 .  38 

15.63 

19  . 

500 

99.60 

2.12         3.95          7.15 

12.29 

1     .-'.3 

.v 

1000 

99.79 

2.17         4.0            7.21 

12.37 

19.73 
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TABLE  9     LONGITUDINAL  HEAT  CONDUCTION  EFFECT 

^in'Snax  M  l-° 

.    Cr/Ctain  =  10.0   ,    (hA)*  =  As* 

«   .5 

•fc 

Conduction  ef£ect;         % 

NTUo 

0 

.015 

.03              .06                  .12 

.24 

.48 

1 

49.96 

0.65 

1.23           2.24               3.84 

6.04 

8.51 

2 

66.59 

0.90 

1.72           3.20               5.67 

9.32 

13.84 

3 

74.91 

1.02 

1.97           3.70               6.62 

11.04 

16.64 

4 

79.89 

1.10 

2.12           4.00                7.21 

12.09 

18.34 

5 

83.22 

1.15 

2.23           4.21                7.61 

12.80 

19.48 

10 

90.78 

1.28 

2.49           4.73               8.59 

14.51 

22.13 

20 

95.11            1.36 

2.65           5.05               9.20 

15.58 

23.76 
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TABLE    10     LONGITUDINAL  I1LAT  CONDUCTION  EFFECT 


Qnln/Cmax  B  LO   >   Ct/^atn  ■  1.0   ,    (hA)*  -  As*  -   .25 


«   % 

Conduction  e 

ffect   £&-  % 

> 

NTUo      A 

0 

.025 

.05              .1 

.2 

.4 

.8 

1 

46.59 

1.02 

1.90           3.38 

5.57 

8.25 

10.84 

2 

60.13 

1.43 

2.69           4.88 

8.22 

12.50 

16.84 

3 

66.64 

1.70 

3.21            5.80 

9.78 

14.88 

20.07 

4 

70.10 

1.92 

3.61            6.50 

10.88 

16.47 

22.16 

5 

73.90 

2  JO 

3.94           7.05 

11.73 

17.65 

23.66 

10 

81.28 

2.74 

5.05           8.81 

14.25 

20.93 

27.63 

20 

86.64 

3.52 

6.27         10.56 

16.50 

23.61 

30.62 
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TABT..E   11      T/>NRITUDTNAL  UEAT   CONDUCTION   EFFECT 


C,ain/Cnax  ■  1.0   ,    Cr/C„,in  -  10.0    ,    (hA)*  =  As*  =    .25 


NTUo 


e    % 


0 


Conduction  effect  ^—  '7C 


.025  .05 


.1 


.2 


.4 


.8 


1 
2 
3 
4 
5 

10 
20 


49.96 
66.60 
74.91 
79.89 
83.22 

90.78 
95.11 


1.07 
1.47 
1.68 
1.80 
1.88 

2.1 


1.99 
2.79 
3.20 
3.44 
3.61 

4.03 


2.23  4.29 


7.92 
5.06 

5.84 
6.32 
6.64 

7.44 
7.92 


5.84 

8.59 

10.00 

10.86 

11.44 


8.68 
13.21 
15.53 
16.93 
17.88 


12.83         20.10 


13.69         21.46 


11.49 
18.05 
21.43 
23.48 
24.86 

28.02 
29.83 
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:  12     L                                             (JCTIOI         '  ECT 
inln/Coax       ■'■.■■     r  La          -      •      ' ■•>"    '   '^  "  '   r 

"  % 

Conviction  effect              '?0 

NTUo 

0 

.01              .02                .04                .08                .16                .32 

1 
2 
3 

« 

in 

20 

vr. 

100 

<*7 .  63 
61.56 
68.46 
72.74 
r.72 

83.20 

88.69 

92.41 

95.50 

0.38           0.73           1.36                2.'.1                 .98 

0.53           1.02           1.93               3.5?             6..                 9.62 

0.62           1.21           2.30               4.24              7.37 

0.70            1.36           2.60                 ..78              3.32            L3.33 

0./6           1.49            ".84               5.22             9.0*                .50 

1. 01            1.96            3.70                 6.70                                  17. 

1.32           2.^2           *.66               3.21           12                 20. 

1.68           3.15           5,^4               9.58           15.32           22.74 

2.27           4.08           6.94             11.17            17.10            24.67 
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•    ■     !  I     I     .  T    HEAT  CO      "•        ON  EFFECT 


t'fUo 


^ia/Cmcx       .9   ,   Cr/C^in  »  2.0   ,    (hA)"*  =  A»*  -  I  .0 


i    % 


Conduction  effect >-      % 


.01 


.02 


.04 


.03 


1 

50.30 

0.40 

0.77 

1.43 

2.53 

2 

66.92 

0.55 

1.07 

2.03 

3.71 

3 

75.30 

0.63 

1.24 

2.37 

4.37 

'♦ 

*QJC 

0.69 

1.35 

2.59 

4.82 

5 

83.85 

0.73 

1.43 

2.76 

5.14 

0 

91.99 

0.85 

1.67 

3.23 

6.04 

.10 


20 


96.82 


0.90 


1.79 


3.49 


6.59 


4.17 
6.41 
7.66 
8.50 
9.10 

10.71 

11.74 


.32 

6.23 
10.16 
12.36 
L3.80 

14.82 

17.44 
19.11 
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TABLE  14     LONGITUDINAL   HEAT  COHDUCTTON  EFFECT 
^in/Coax  -   .9    ,   C^in  =  10.0   ,    (hA)*  =  As*  =  1.0 

*% 

Conduction  effect,  —  % 

NTTIo 

0 

.01              .02             .04               .08                  .16               .32 

I 
2 
3 
4 
5 

10 

20 

51.22 
68.81 
77.67 
82.98 
86.52 

94.36 

98.36 

0.41           0.78           1.45             2.56               4.23             6.32 
0.56           1.09           2.06             3.77                6.49           10.28 
0.63           1.24           2.38             4.39               7.68           12.38 
0.68           1.33           2.56             4.76               8.39           13.65 
0.70           1.39           2.68             5.00                8.87            14.49 

0.75           1.48           2.88             5.48               9.87           16.35 

0.65           1.33           2.70              5.35             10.03           17.08 
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TABLE   15     LONGITUDINAL 

HEAT  CONDUCTION  EFFECT 

Cmin/Ctu 

ox  ~   *9   , 

Cr/^min 

•1.0,    (hA)*  -  As*  = 

.5 

! 

e   <V 

Conduction  effect              <> 

\ 

NTUo 

X 

0 

.015 

.03 

.06                  .12 

.24 

.48 

1 

47.53 

0.54 

1.03 

1.88               3.23 

5.U 

7.23 

2 

61.40 

0.74 

1.43 

2.67               4.76 

7.89 

11.81 

3 

68.26 

0.87 

1.69 

3.17               5.70 

9.55 

14.48 

4 

72.51 

0.97 

1.88 

3.55                6.40 

10.73 

16.29 

5 

75.47 

1.06 

2.05 

3.86                C.94 

11.63 

17.64 

10 

82.98 

1.35 

2.61 

4.89                8.69 

14.33 

21.38 

20 

88.34 

1.68 

3.23 

5.97             10.39 

16.70 

24.35 
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"  "      (  CTUDINAL    '. :.  '  nXJCTION  EFFECT 


Cmin^x    :   -q    »    Cr/Cmin  *  2.0  ,    (hA)*  -  As*  =   .5 


e  % 

Conduction 

effect    Af 

> 

% 

NTJJo 

o 

.03 

.12 

.21 

.30 

1 

50.27     • 

L.oa 

3.40 

4.96 

c.09 

!\ 

66.88 

1.53 

5.08 

7.71 

9.73 

3 

75.23 

1.77 

6.01 

9.23 

11.76 

4 

80.32 

1.93 

6.63 

10.24 

13.08 

5 

83.77 

2.05 

7.09 

10.97 

14.03 

6 

86.28 

2 .  14 

7.45 

11.53 

14.75 

7 

88.18 

2.22 

7.74 

11.93 

15.32 

3 

89.  6  9 

2.29 

7.97 

12.34 

15.78 

9 

90.90 

2.34 

8.17 

12.64 

16.17 

10 

91.91 

2.38 

3.34 

12.90 

16.49 

12 

93.^7 

2.45 

8.60 

13.31 

17.00 

14 

94.62 

2.50 

3.80 

13.61 

17.39 

16 

95.50 

2.54 

8.96 

13.85 

17.67 

18 

96.20 

2.56 

9.07 

14.03 

17.90 

20 

96.76 

2.57 

9.16 

14.18 

18.09 
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!  VL    IFA'E   COIID  ••      ... 


■     0  ;,     .      .    in  o       /*,»>*  a  An*  ■     5 


e  % 

Conduction  effect 

7* 

\ 

NTUo 

0 

0.015 

o.o;. 

0.06 

0.12 

0.24 

3.48 

1 

51.22 

0.58 

1. 10 

2.01 

3.45 

5.45 

7.73 

2 

68.80 

0.80 

1.55 

2 .  89 

5.15 

8.55 

12.83 

3 

77.66 

0.92 

1.78 

3.35 

6.04 

10.10 

15.  5* 

4 

82.98 

0.98 

1.91 

3.62 

6.58 

11*18 

17.22 

5 

86.51 

1.02 

1.99 

3  79 

6.39 

11.85 

18.33 

10 

94 .  36 

1.07 

2.12 

4.11 

7.65 

13.29 

20.77 

20 

98.35 

0.96 

1.96 

3.97 

7.71 

13.80 

21.91 
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TABLE  18     LONG        D   HAL 

HEAT  c          TIC  :  :  FFECT 

Stin  W 

■       Q 
X           "      * 

Cr/CB»i  , 

1.0  ,    (hA)*       Kb* 

.25 

e     % 

Conduction  effect              <? 

-TUo 

0 

0.025 

0.05 

0.1                  0.2 

0.4 

0.3 

1 

47.60 

0.37 

1  .  62 

2.90               4.82 

7.21 

9.57 

2 

61.29 

1.18 

2.26 

4.13                 7.08 

10.98 

15.09 

3 

68.13 

1.39 

2.65 

4 . 87                c .  39 

13.10 

13.12 

4 

72.37 

1.54 

2.94 

5.41                9.32 

14.54 

20.11 

5 

75.31 

1.66 

3.17 

5.83              10.03 

15.61 

?.l .  5A 

10 

N2.74 

2.05 

3.92 

7.17              12.16 

18.61 

?5.35 

20 

33.02 

L 

2.43 

4.67 

G.45              14.04 

21.04 

2J  .16 
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. 

B  1                         ...  ^L  HKM      .-   Dt    riON 

CI 

1 

,   G   '■     in   -   LO.      ,    .    .  *       Kb* 

-   .25 

c   % 

Conduction  effect,          % 

0 

\/ .  .  ^  .- 

0.05            0.1                   0.2 

0.4 

0.8 

- 

J  J 

3 

1.74            J.  10                >.l  > 

. 

2 

6*  ,  80 

1.29 

2.^5           4.49                7.72 

K.Ct 

16.7  ■• 

3 

7  7.'' 

1.47 

2.82              .21                 \06 

14.32 

20.12 

4 

82.98 

1.58 

3.03           5.64                   37 

1  5 .  7C 

22.  V 

5 

86.51 

1.64 

3.17            5.93             10.41 

16.62 

23.5 

10 

.3' 

1.74 

3^1            6.    8              11.56 

18.67 

26.62 

20 

98.35 

1.59 

i.26           6.'  6             11.89 

19.5/ 

28.  OC 
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TABLE  20  LONGITUDINAL  HEAT  CONDUCT  J.  ~l  FECT 


Cain/Cmax  ■  0.5  ,  Cr/Cctn  -  1.0  ,  (hA)*  *  As*  -  1.0 


«    fc 

Conduction  effect 

)           /o 

NTUo^ 

0 

0.01 

0.02           0.04 

O.Cr 

0.16 

0.32 

2 

66.89 

0.35 

0.68           1.29 

2.37 

4.15 

6.71 

3 

74.01 

0.41 

0.80           1.53 

2.84 

5.02 

8.25 

4 

78.16 

0.45 

0.89            1.71 

3.19 

5.65 

9.29 

5 

80.89 

0.49 

0.97           1.86 

3.46 

6.12 

10.06 

10 

87.15 

0.05 

1.26           2.39 

4.38 

7.57 

12.15 

20 

01.07 

0.86 

1.64            3.04 

5.37 

8.94 

13.87 
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.     1  :•   21     l             Oil  ...    HEAT  C0NDUC1  r.;  1  EFFECT 

' 

Qailn'^nax  "  ^-  5  » 

Cr/Cttln  -  10- 0  ,    0*)*  "  /s' 

=  1.0 

e  7o 

Conduction  effect ;           ?o 

NTUo 

I    ° 

0.01 

0.02            0.04              0.08 

0  3  I 

0.32 

I 

56.42 

0.28 

0.53           0.98              1.74 

2.90 

4.37 

2 

77.  34 

0.37 

0.71            1.35              2.49 

4.36 

7.06 

3 

87.30 

0.38 

0.74           1.43              2.68 

4.82 

8.03 

4 

92.60 

0.35 

0.70            1.36              2.60 

4.80 

8.29 

5 

95.59 

0.21 

0.62            1.22              2.39 

4 .  sr- 

8.1 

10 

99.62 

0.1.0 

0.22              .49              1.14 

2.68 

5.96 

20 

99.99 

0.01 

0.01           0.05             0.23 

1.02 

3.65 

41 


i 

lain'   tiax       lJ,J   » 

Cr/Crain 

-  1.0   , 

(hA)      -  Aft 

-0.5 

-  7o 

Conduction  effect^-   -% 

m  *  x\^ 

0 

o.oi 

0.02 

Q.04 

o.o; 

0.16 

0.32 

i 

50.97 

0,2? 

.53 

0.98 

1.71 

2.75 

3.99 

2 

65.65 

0.35 

.69 

1.32 

2 .  40 

4.15 

6.56 

3 

?2.4A 

0.30 

.78 

1.49 

2.78 

4.92 

8.0A 

4 

76. A3 

Q.A2 

.87 

1.62 

3.04 

5.46 

9.0A 

5 

79.10 

0.45 

0.89 

1.72 

3.ie- 

5.86 

9.79 

10 

8S.'+4 

0.56 

1.09 

2.12 

4.00 

1     8 

11.95 

20 

89.72 

0.74 

l.AA 

2.  72 

* 

A. 95 

8.58 

13.90 

42 


\ 

•     i  ' 

•    .>       . 

/, 

n  n 

! 

(hA) 

Q 

•X 

Conduction  effc 

t 

9 

,.-  to 

0J 

5.02 

0,0& 

o.^a 

lb 

0.32 

— 

1 

1 

v.    M 

! 

61 

L.?l 

I  .• 

4.54 

77.32 

0.4 

0.83 

i  »:7 

2 

4 

7.7^ 

'! 

87.27 

0,44 

0.37 

1 .  68 

3.1^ 

■    61 

9 .  i  2 

4 

37.57 

0.42 

0.82 

1.62 

3.10 

|      v      « 

9.73 

95.56 

0,3? 

0.74 

1,47 

2,90 

5 . r.  3 

o    /*. 

10 

99 .  61 

o.>^ 

0.22 

0.63 

1  .54 

3.^ 

8.13 

2< 

9^   99 

0..  ■ 

0.02 

0.08 

).43 

1.92 

'..nt. 

1ABLE   24 

LONGITUDINAL 

C  h  DrCI-r.N   ti-Tt 

XT 

^nin'Sa 

■  0 

tax 

.5 

*  Cr'^iBin 

-   1.0    , 

(HA)*  =  As* 

-  0.25 

•  ?o 

Conduction  effect  ,  =■ —  % 

> 

NTUo 

0 

.01 

.02 

.04 

.  )8 

.16 

.32 

1 

50.21 

o.j/ 

0.65 

1.13 

2.01 

3.12 

4.:4 

2 

64.39 

0,3;} 

0.75 

1.42 

2.61 

4.47 

6.89 

3 

70.90 

0.39 

0.77 

1.51 

2.89 

5.19 

6.38 

4 

74.91 

0,39 

0.79 

1.58 

3.09 

5.70 

9.43 

5 

77.57 

0.40 

0.81 

1.64 

3.26 

6.11 

10.24 

10 

04.12 

0.49 

0.99 

1.99 

3.97 

7.53 

12.71 

20 

88.72 

0.66 

1.3^ 

2,54 

4.92 

9.07 

14.97 
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TABLE  25     LONGITUDINAL  HEAI   CONDUCTION   EFFECT 
<W<W  *  0-5   ,   Cr/CWn  -  10.0    ,    <hA>*  =  As*  -  0.25 


e  7c 


Conduction  effect 


>  " 


NTUo 


.01 


1 
2 
3 
4 
5 

10 

20 


56.40  0.43 
77.30  0.58 
87.24  i  0.60 


92.53 


0.57 


95.53   0,50 


99.60   0.18 


99.99  ,  0.01 


.02 


0.41 


0.04 


.04 


1.01 


0.22 


.08 


2.59 


1.13 


'.6 


6.18 


4.29 


.32 


0.81 

1.47 

2,51 

3.92 

5.46 

1.11 

2.10 

3.79 

6.36 

9.59 

1.18 

2.27 

4,24 

7.42 

11.66 

1.12 

2.22 

4.28 

7.77 

12.62 

1.02 

2.05 

4.10 

7.73 

12.99 

12.19 


10.59 
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< 


CD 


5       6      7      8       0.01 
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Vc 

mi 
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M 

1 

,^ 
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% 

^ 

^x 

X" 

^ 

20 

0 

> 

3 

2 

10 

. 

tlCr/rmin=2 

^ 

1 

0 

NTUo 

1 

10,20  ! 

3 

i 

in 

1 
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APPENDIX  1. 

Derivation  of  iterative  equations. - 

Equattons  (1)  and  (2)  from  page   can  be  put  in  the  following  form: 


*  (, '-,  l )   +  *  ■  • 


'Xti,j) 


(1.1) 


—  T  fT7      —  T       C 


'X{.i,fi  ;/  +  '■ 


(1.2) 


where 


= 


and 


- 


. 


Equations  (1.1)  and  (1.2)  may  be  written  as 


- 


' 


(1.3) 


and 


' 


' 


(1.4) 


- 


■') 
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uy   ;,.U!ing  to  equation  (1.3)  the  product  of  (1+C.)(1.4)  #nd 


I 


(4) 


where 


and  i__ 


From  equation  (1.1) 


• 


(6) 


where 


= 


■ 


• 


am; 


For  the  elements  from  the  first  and  last  row  equation  (1.2)  becomes 


(1.5) 


an<i 


• 


(1.6) 
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respectively.   I,  to  same  procedure  as  before  and  solving  for  the  matrix 
outlet  temperatures,  the  following  equations  are  obtained: 


'  ' 


and 


(3) 


.. 


(5) 


for  the  first  and  last  row  elements  respectively,  where 


, 


and 


D.- 


Equations  (7)  and  (8)  may  be  written  as 


r 


- 


(1.7) 


and 


— 


(1.8) 


Equation  (1.8)  reduces  to 
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and 


<-p     t    t1 
■  |)  +  W^.)""  *<*-'»j)~  W      (1.9) 


_  p  fq1 


for  the  first  and  last  vow  elements  respectively,  where 


(1.10) 


=  N  I  L 


I 


By  comparison  of  equations  (1.7)  and  (1.8)  with  equations  (1.1)  and 
(1.2)  it  follows  that  the  form  is  the  same  and  only  the  coefficients  are 
different.   The  same  holds  true  for  equations  (1.9)  and  (1.10)  with 
equations  (1.6)  and  (1.5)  respectively.   By  substitution  the  following 
equations  are  obtained: 


_  r 

.   ■   "   .     .  ,)  l 


T        T 


(9) 


. 


' 


(10) 


~ 


(11) 


for    the    first,   middle  and  last  row  elements  respectively  and 
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(12) 


which  is  the  same  for  all  the  elements,  where 


> 


and 


i 


5A 


Con     -  .  •       i.- 

Kqu.ttion  (1 .2)   car»   o-?  y\SL   In   *'..<     following    roria: 


<         - 
n   the  energy  balance  equations,    the  aide  ot  Crnax  was  assunv  i  to  be  the 

hot  sitfc   so  that  by  the  second  law  of  thermodynamics 


and    from  equation   (1.11) 


■ 


■ 


Reari&n;  ing   equation  (1.12) 


(1.12) 


, 


(1.13) 


Again  by  the  second  law  oi  thermodynamics  the  matrix  temperature  variation 
in  the  direction  of  fluid  flow  is  such  that 


(1.14) 
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..vij.i 

(1.15) 
so  that 


(1.16) 
Substituting  equation  (1.16)  in  equation  (1.13) 


(1.17) 


Three  cages  are  investigated: 

1.-   If  C.  C,  then  equation  (1.17)  becomes 

2.-   If  C^  *  0  then  equation  (1.17)  becomes 

0  ■ 

3.-   If  C.  ■  2C2  ■  0  then  equation  (1.17)  becomes 


From  the  above  casr.s  it  can  be  seen  that  a  sufficient  condition  for 
equation  (2)  to  converge  is  that 


Substituting  the  values  for  the  constants  and  rearranging  yields 
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For  the  first  and   U;St  row  elements  the  conve  gence  criterion 
becomes 


which  for  substitution  of  the  parameters  gives 


— — 


This  condition  ts  loss  stringent  than  for  the  middle  row  elements, 

Similarly,  for  the  side  o   Coin  the  sufficient  condition 
for  the  middle  row  elements  becomes 


which  for  substitution  of  the  parameters  gives 

-   . 

1 

For  the  first  and  last  row  elements  it  is  found  that 


which  requires  that 


— 


Is  the  sufficient  condition.   This  again  is  a  less  stringent  condition 
than  that  imposed  by  the  middle  row  elements. 
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Differentia  s.  - 

Hausen       4       has  fclv.n  a   theory    fur  rcgtn retort  without  heat 
conduction.      .he*    the  bat  t  uss-tf  t  «  .&  ,i>  i.^ade,   but  tb»    Longitudinal 
conduction  of  heat   is  included,    (See  page       )    the  differential 
equations  obtained  are 
for  the  side  of  Cn&x  (gas   flowing   frcro     ;    ■»  0  to         ■  L) ; 


- 


and  for  the  side  of  Cmtn   'reversed  gas  flow) 

-  - 

A  temperature  scale  can  be  us.d  for  *  hich  the  gas  entrance 

temperature  at  the  side  ol  Cri^n  is  zero  and  that  at  the  side  of  C 

pax 

is  unity.     The  boundary  coi.a.tions  are  then 


and 





APPENDIX  2 

Computer  program. - 

The  Fortran  system  was  used  to  program  the  problem  for  the  CDC  I60'f 
digital  computer. 

The  sir.   diraensionless  parameters  are  used  as  input  together  with 
increments  and  factors  to  change  the  values  of  NTUo,    and  As,  so  that 
solutions  can  be  generated  for  several  combinations  of  the  parameters 
without  having  to  feed  into  the  computer  these  different  sets  of  para- 
meters every  time. 

Another  input  is  the  number  of  subdivisions  of  the  three  streans 
together  with  the  increments,  so  that  when  the  computer  finishes  the 
first  series  of  runs  for  a  given  number  of  subdivisions  it  automatically 
increases  the  subdivisions  by  the  specified  increments  and  the  two  values 
of  effectiveness  obtained  are  extrapolated  to  an  infinite  number  of 
subdivisions. 

In  order  to  obtain  the  initial  estimates  of  matrix  temperatures, 
the  problem  is  first  solved  neglecting  the  longitudinal  heat  conduction 
as  explained  previously.   For  this  case  the  equations  derived  in  Appendix  1 
are  simplified  to  (for  detailed  derivation  see  Reference  I) 
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whtre 


I    + 


4- 


C, C, 


-I 


= 


Nr 


-I 
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The  inlet  matrix  temperatures     )  necessary  to  Initiate  the 
problem  for  zero  heat  conduction  are  estimated  by  assuming        , 
(See  Fig.  lb)  md  feeding  this  value  into  the  computer.   The  remaining 
inlet  matrix  temperatures  are  obtained  from  the  relation 


Mr 


The  number  of  passes  necessary  to  meet  the  desired  accuracy  in  the 
solution  changes  accordingly  with  the  initial  estimates  of  matrix  inlet 
temperatures,  but  considering  that  the  computer  time  required  to  work  the 
problem  for  zero  heat  conduction  is  so  small  as  compared  to  the  case  with 
conduction,  that  these  estimates  do  not  become  important.   Only  for  high 
values  ot  Cx'/Cmin  these  estimates  may  be  of  significance  since  for  this 
case  the  change  in  the  reversal  condition  after  every  pass  is  very  small. 

In  general,  it  was  found  that  an  estimate  T.  ,^  jv 
of  0.4  to  0.6      for  ^/(^^^S.O 

and         0.7  to  0.9      for  Cr/C„,in^5.0 

together  with  the  above  relation  for  Tx , <  ,v  should  converge  to  the 
solution  in  a  reasonable  time. 

The  temperature  distribution  is  calculated  by  repetitive  use  of  the 
iteration  equations. 

The  effectiveness  may  be  defined  as 


[Wi,    "     '-Xp/  Lrio  ~^'r><- 


C"l*~      {t%i—t%JL)  ^**  I-Tvi, 
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ol 


M«   A  ^ 


^x»     ~n7  2.    ■- 


0 


It,,     ^  /_  Tu(   .,  -  U 

■ 
therefore,  the  expression  for  effectiveness  for  the  definition  reduces  to 

p  =  —  y  t 

which  is  the  expression  used  for  the  computation  of  the  effectiveness. 
The.  heat  balance  error  is  computed  from  the  relation 


error  =  ■ — 


which  reduces  to 


error  = 


- 
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The  e.:trapo1  ".ess  i-  com  the  relation 

(v.     •  "i'-.  -> 

e ... 


where 


RE  =    •' 


. 


The   input  parameters  in  Fortran  notation  are 

Pl   =  C,ain/Cmax   >   *2  "  V^in   •    1>3  "   <llA>*   » 

p4  =  WTUo    ,    P5  -    \    ,    ?6  =  As* 

The  initial  values  of  NTUo,  \  and  As*  are  changed  by  the  following 
relations 

PA  =  p4  -:-  P4IN 

P5  «  P5*P5FA  +  P5IN 

P6  »  P6  -  P6IN 
where  P4XN,  P5FA,  P5IN,  P6IN  are  also  input  values.  The  inputs  P4FI , 
P5FI  and  P6FI  are  the  final  values  of  P4,  P5,  and  P6  respectively. 
The  program  is  written  so  that  the  computer  after  working  the 
initial  set  of  parameters,  i'  re  duces  ~G  by  P6IN  every  time  until 
P6  becomes  less  than  P6FI .   Then  it  sets  P6  equal  to  the  initial  value 
and  changes  P5.  The  procedure  is  again  repeated  until  P5  becomes  greater 
than  ?5F1.  Then  P5  is  set  equal  to  the  initial  value  and  P4  is  increased 
by  P4IN.  The  cycle  is  repeated  until  P4  becomes  equal  to  or  greater  than 
P4FI. 
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■ 
.  i  Lecisions 

"8  Or  '        '  ■  ■  '  )  '  I •:  ' 

a     and  7b  sho*>  th<  comi      Low 
i  wl        ;-    I  by  the  For  tan  listing  of 
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1      INPUT   VALUES 


1 l , 

I  2      COMPUTE   INITIAL  Tx(i   X)  1 

|  1 ^ 

13      SET   tx(lij)    =   l.o     AND      tn(l 


M  =    PASS   Nj.--.8ES 

IC   =   1,    INITIAL  SUBDIVISIONS 

10   =   2,    INC.1E  MJTED   SUBDIVISIONS 


)    =    0 


NO  CONDUCTION  CASE 

COMPUTE  Tx(l,j)  ,  TD(r,g)  ,  e  *"D  ERROR 


CONDUCTION  CASE 

23 

COMPUTE  CONSTANTS 

2» 

SET  M  =  1 

25 

ITERATION  ON  SIDE  OF  CMX 

26 

ITERATION  ON  SIDE  Or  C„in 

2? 

COMPUTE  e  AND  EBROS 

FIO.  7a     OENEKAL  FLO*  DIAGRAM 
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COMPUTE  CONSTANTS 

5 

SET  M  ■  1 

t 

'  c„ax 

7 

ITERATION  ON  SI  :>E 

8 

COKPUTE  e  »ND  ERROR 

11      STORE 

Tx(i,l),e,    ER-fOR 


13      STORE  Tx(i,j),    Tn(f,g) 


H18      REVERSAL  CONDITION 
BY   T"Cr.N„*l) 


,32      INCREASE   Nr  ,    Nx,    N„     I (T) 


riG.    7o  f'UMl    DIA,.<AM   FOR   HO   CONDUCTION 
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PROGRAM 

REMARKS 

PROGRAM    BAHNKE 

DIMENSION    THlUOiMOltUHIUO.UOl ,TC(«O.UOI ,UC(UO,UO) , ERR  1921  ,  TCOI92), 
lUMCCt0,UO),UCC(tO,i»0l  ,UH  1  IUO)  ,UH2(U0)  ,EFI20,20),FR!20.20) 

1  F0RMATI6FS.2I 

READ     l,PI,P2,P3,Pll,PUIN,  PUf  1 

2  F0RMATI6I5I 

READ    ?,NM,NH,NC.NMI,NHI,NCI 

III  POT 

NM   •    Mr.    NB    ■    Ni ,    NC    .    Is 

J            FORMATIFS.3I 

READ     S.UMI 1,1) 

EM*NM 

UO     UOI     K  = 1 ,NB 

t=K 

MOI       UHIK, 1 I=CIEM+1,0-EI»OHI1 , 1 1 1/EM 

INPUT!   Did ,19   -  ta(x,i) 

ESTIMATES   OF   Tx(i,j) 

4            FORMAT17F7.2I 

READ    It.PS.PSIN.PSFA.PSf  1  ,  P6,  P<>  1  N  ,  P6f  I 

INPUT 

SS'PS 
S6=P6 

STORE    INITIAL   PS   AND   P6 

5            NM2=NM»NM! 

NH2=NH«NHI 
NC2=NC«NCI 

NUMBER   Or   SUBDIVISIONS   TOR    2nd    RUN 

NE  l=NM«INH»NC 1 
RE  l=NE  1 

NE2=NM2«(NH2*NC2I 
RE2=NE2 
RF=RFI/IRF?-REII 

rfE    •    VARIABLE   USED  TO   EXTRAPOLATE   EWECTIVENES* 

EM2=NM2 
U06    L=  I.NM2 
t=L 

RT»1.0-CIE-I.0I/EM21 
6            UH2IIUUHI  1,  ll«RT 

ESTIMATES   OF  T»(l,l)    FOB    2nd   RUN 

7            FORMAT!  //HtHNO    OF    ELt  MEM  r  S3X,  31  S  ,  l«X  ,  31  S  ) 

WRITE    OUTPUT    TAPE    2 , 7 ,NM ,NH, NC , NM2 , NH2 ,NC2 
PRINT    7,NM,NH,NL,NM2,N.H2 , NC 2 

8           IC=  1 

INDICATES    1st    RON,    USING    INITIAL   No   OF   3UBDIV. 

9            EM=NM 
EH=NH 
EC-NC 

rlXED   TO    FLOATING    POINT   CONVERSION 

NMH-NM- 1 
NHPiNHt 1 
NCP^NC* 1 
NMPiNM* 1 

DEFINITION    OF   NE*    VARIABLES 

DO     10    J=I,NH 

10  THCI,J!=1.0 
UO     11    L-l.NC 

11  IC 1  1  ,L 1=0.0 

SET   'x(i.j)    ■    1-0      AND      tn(l,g)    -    0.0 

12          AlO=l.0*EH/(PI«^2«EHI»(2.0«6M»P3)/IPU#Pl»(l.0»P3)l 
AI=2.0/AIO 

A2D«l.0»(Pl»P2»tHI/EMt(2.0«EH»P3«P2l/(PU»(l.OtP3)l 
A2*2.0/A20 

A40*I.Q«EM/IP2«rC)     H2.0«EMI/(P"l»11.0»P3ll 
A3»2.0/A3D 

AMD=I.0*IP2«ECI/EM    ♦(2.0«EC»P2)/(PM«(I.0+P3I> 
A"l  =  2.0/AU0 

CONSTANTS   roR    THE    NO   CONDUCTION    CASE 

C               HOT    SIDE    NO    CONDUCTION 

121         M*l 

M   DENOTES    PASS    NUMBER 

13         00    IS    J«ltNH 

DO     IS    1*1, NM 

THU«I,J»  =  TH(!.J)»AIMUH(I,J>-TH(I,JI) 
IS          UHll.J* 1 l=UH(I , J|-A2«IUH( 1 , JI-THI I , Jl) 

ITERATION    ON    SIDE   Or    Caa, 

C               COID    SIDE    NO    CONDUCTION 

00     17    K-l,    NM 
17         UCIK , 1 1*UH(NM» l-K.NHt 1 1 

DO    19    l=l,NC 

0019    K=  l,NM 

TC(K»  l,LI*TC(K,L  l»A3«(UC  (K.LI-T1.  IK,  I  I  1 
19        UCIK.L*  M=UCIK,U-AM«(UC  (k,l>-TC(k,l>) 

ITERATION   ON    SIDE   Or    C.ln 

GUO.O 

DO    20    J=1,NH 
20         GIMH(NM»I,J l«GI 
TH0=GI/tH 

TBO   •     Lt«(Nr.l,J)j    srg 

sr 


on    2  1    L=I,NC 

G2=TC(NM*1,L HG2 
TCO(MI=G2/EC 


ERR ( M ) = I .0- ( P I • TCO ( M I  1 / I  I  . 
IFIERRIMI     123,2li,?li 

23  ERRIM)*-ERRIMI 

2"<  IH  fcRR(MI-l|.0f-6)     30,30,25 


TCO   =    EFFECTIVENESS 


25 
26 

21 
28 
29 


IFIVO-M)    26,26,215 


COMPUTE   HEAT   BALANCE   ERROR   AND  COMPARE    II 
HITB   MAXIMUM   SPECIFIED 


If ( tRRlM-1 l-tRR(H) )    2V.2V.27 

PAUSE     1  I 

IFISENSfc    SWITCH     1 )     32,  121 

PAUSE  12 

00  TO  28 


CHECK  CONVERGENCE 

PAUSE  11  =  SOLUTION  CONVERSES 

PAUSE  12  =  SOLUTION  DIVERGES 


215   DO  22  1=1, NM 
22    UH( I, I )=UCiNM+l-I,NC*l ) 
M=M»1 

GO  TO  13 


REVERSAL  CONDITION 


SOI 
302 


T2=TC0(M-I l-TCO(M) 

IFIT2)     301,302,302 

T2=-T2 

If  (T2-''.0E-6I     305,305.25 


CHECK    DIFFERENCE    IN    EFFECTIVENESS    IN 
THE   LAST  T*C   PASSES 


305 

51 

312 


IF(  IC-1)     31,31, 33 
L)0     312    J=1,NM 
UH1I J)=UH( J,  I  I 
EFI  I  ,  1  )  =  TCOIM) 
ER(  1,  I  )  =  ERR(M) 
GO    TO    3V 


STORE  Txli.,1)    •    EFFECTIVENESS   AND 
ERROR   FOR   THE   1st   RUN 


F0RMAIIUF8.2, I5,F10.ft,FV.?,I3) 

PRINT    32,P1,P2,P3,PU,M,TC0(M).ERR(MI  ,IC 

PAUSE     13 

IFISENSE    SWITCH     1 )     321  ,  l?l 

1FISENSE    SWITCH    2)     1,322 

NM=NM2 
NH=NH2 
NC=NC2 
GO     TO    '.  


EFF=TCQ(M)«1 1.0*RE)-EF( 1 , I )«Rf 

Hl  =  tFF _  

FORMAT! //UF8. 2. l6X,2(FIU.6,r«.2l,F10.6) 

WRIIE  OUTPUT  TAPE  2,3>I,P  1  ,P2  .P3.PU  .EF(  1  ,  1 1  .ER(  I  •  1  I  .  TCO(M)  ,EKR(M| 


INCREASE  NO.  Of  SUBDIVISIONS  AND  REPEAT 
1st  RUN 


EXTRAPOLATE  EFFECTIVENESS  AND  STORE 


35 
39 


35  L=1.NM 
UH2IL )=UH(L,  I  I 


STORE  Txd,!)  FROM  2nd  RUN 


IA  =  2 
IB=I 


HO 
U  I 
U2 
U3 


DO    UU    1=1, NM 
00    U2    J= l.NHP 
UHC| I, JI=UHI I , J) 
DO    UU    1=1, NOP 
UCCI  l,L)=UC(  1,1  I 


SUBSCRIPTS   FOR   ERROR   AND   EFFECTIVENESS 


UHC(I.J)    =   Tx(lij)    ,    UCC(I.L)    =   T0(fig) 
FOR   CONDUCTION   CASE 


CONSTANTS 

E1  =  IPH«(  1.0*P3)  1/(2.  Clr  I 

El  1=  I.0*E1 

T2=IP2«ECI/EM 

E3=0.5«P5«EM 

EU=E1+(E  11*12. 0>E3*E2> ) 

E5»EI+(E11*(E2+E3I ) 

F1=(2.0*El)/f4 

F2=(EUIElli|2.0.fi-f?ll  | /EM 

F3»(E3*E 1 1 l/EM 

FM=(  1.0-E1 l/E  1  I 

F5=EI/E1 1 

F6=I2.0*E1 l/t5 

Ff-lEI+lEl 1*(E3-E2)))/E5 

F8=(E3*E 11 r/ES 

CI=(P1*Pli.(|.0*P3l l/(2.0*P3*EM) 

CI  1=  I.0»C1 

C2=1PI*P2»EH)/E~ 

C3=(P1*P5«EM)/(2.0*P6) 

CU>CH(Cll*(C2*2.0*t 

c5=cmc i i*(C2*c  ii  i 


CONSTANTS   FOR   CONDIJCTKN    CASE 


69 


D  I  =  t  2 . 0«C I 1/Ot 

02=  (C  t  ♦  I  C  1  1*1  2.0«C  5-C?  I  I  I/O. 

03=1C3«C  I  I  I/O* 

DU-1  1.0-C I l/C I  I 

05=CI/CI 1 

D6=(2.0«CI)/C5 

07=(CU(CI  I»IC3-C2  I  I  I/C5 

D8=(C3«C1 I I/C5 


CONSTANTS  FOR  CONDUCTION  CASE 


5  3    M=t 

C      HOI  SIDE  WITH  CONDUCTION 

60    00  6?  JM.NH 
1=  I 

6  1    UHCI  l,J«ll=U6«TH(  I, JI-D7«UHC(  1 , J  I »D4» (UHC12, J)»UHC(2,Jtlll 

TH(2,JI  =  Dl4«TH(  l,J)+05«(UHCI  I,  JI»UHCI  I  ,  J*  I  I  I 
62    U06U  1  =  2.  NMM 

UHCI  l,Jtll  =  UUTHII,  J  I  - 02 • UHCI  I  ,  J  1*04  •(  UHCI  I-  I  ,  JI+UHCI  1-1  ,  Jt  I  I* 
IUHCI  !«  I.JI+UHCI  I*  I,  J*  I  I  I 
6»    THll»t,Jl  =  DU«THII,J)«D5»ll>HCII,J)tUHCII,J»ll) 

UHCINM, J« 1I«U6«TH(NM,JI-C7»LHCINK,JI »D8« 1UHC ( NMM , J  I ♦UHCINKM, J* I  I  I 

66  THINMP,  J)=£)li«THINM,  J  H-D5  •  I  UHC  I  NM  ,  J  It  UHCINM,  J*  I  I  I 
L  =  L»  I 

IFI3-LI  67,67,61 

67  CONTINUE 

68  DO  69  K=I,NM 

69  UCCIK.  I  )=UHCINMP-K ,NHP I 

L      COLD  SIDE  WITH  C0N0UCTION 

70  DO  7b  L=I,NC 
J=l 

n         UCCl  1iL+I)»F6»TC(  ltL)-F7»UCC(  1 , L I *Ffl» ( UCC I  2 , L l+UCC I  2 , L* I ) I 

TC(2,L)  =  FI|«TCI  l,L)+F5«(UCCI   I, LI  ♦UCCl  1  ,L»  III 
72  00    73    K=2,NMM 

UCC(K,L+lt=Fl»TCIK,L>-F2«UCCIK,L)+F3*IUCC(K-l,L)*UCC  (  K—  1  »  L  »  I  I  » 
I UCCIK* l,L l*UCC(K*l,L» I  I  I 
/3  TC(K»1,LI  =  FMTCIK,L>*F5.(UCUK,L)*UCCIK,L*III 

UCCINM.L*  I  l=F6«TC(NM,L)-F7.UCC(NM,LI «F 8» I UCC I NMM ,L I »UCC I NMM , L ♦ I )  I 
7U  TCINMP,L1*FU«TCINW,LH-FS«IUCCINM,L >♦ UCC I NM, L • I  I  I 

J=J*  1 

IFI3-JI    75,75,71 

75  CONTINUE 
(1=0.0 

DO    76    J=I,NH 

76  GUIHINMP,  J  l»GI 

THO=GI/tH 

G2*0.0 

00    78    L=I,NC 

78    G2=TC(NMP,L>»G2 
TC0[MI=G2/FC 


BO 
81 


ITERATION   OJi   SIDE   OF   C«,x 
TWO   PASSES   PER   COLUMN 


ITERATION   ON   SIDE   Or   Cain 
TWO   PASSES   PEB   COLUMN 


ERR(M>=  1.0-IPUTCOIM)  )/(  I.O-THO) 
IFIERRIM))    80,80,81 
ERRIM!=-ERR|MI 
|FIERR(Mt-l|.0E-6)    90,90,62 


32 
83 


BUI 
8U2 
85 


IFI90-MI    83,83,85 

IFIERRIM- I 1-tRRIMI I     HU,8"i,87 

PAUSE    22 

IFISENSE  SWITCH  II  Bi)2,53 

IFISENSE  SWITCH  21  BH,88^ 

M=Mt  I 

00  86  1=1, NM 

UHCI  I  ,  U=UCCINMP-I,NCP) 

GO  10  60 


PAUSt  21 

GO  TO  8UI 

IFIIC-II  681,881,882 

EFI  IA, lBt=0.O 

ERI IA, 181=0.0 

GO  TO  97 


FORMATI  I IHSOL  DIV  FOR  1 3,  I 8X , £9.2 , F8. 2 , 2  I F I  0. 6 , F 9  .2  I  I 

WRITE  OUTPUT  TAPE  2 ,  882,  IC ,  M5,  P6  ,  FF(  I  A ,  I  B )  ,  ER  I  1 A  ,  I II I  ,  TCO I  M)  ,  r:RR{  M I 

go  to  9-Z . 


b^ 


885 
89 


F0RMAT(l|F8.2,E8.2,f  H.2.2X,  I  3,F  10. 6, E 9. 2) 

PRINT  BH<l>PI,P2,P3lP>tfP5,P6,lC,Tr.0(M),FRR(M> 

PAUSE  23 

IFISENSE  SWITCH  1)  885,  ">  3 

IFISENSE  SWITCH  21  I  14,86 

1F(  1C-I)  891, 891. V3 


91 
V2 


'V7 
971 
972 
973 
9711 
T02 
131 


EF( lA.IBI'TCOlM) 
ERI  IA,  IP.)  FKKIH  I 
GO    TO     102 


TWTCOIM-I  I-TC01M) 
IFIT1)    91,92,92 
Tl=-T 1 
IFI  T  I-U.0C-6I     89,89,82 


EFF=TCO|M>.( I.Q.RE l-FFI 10, If  I .RE 
B2M Iflt-EFFI* I00.OI/R 1 


FORMAT  I  J2X,e8.2,FB.?,2(F10.(.,F9.2),FI0.6,F8.3l 

WRJTE  OUTPUT  TAPE  2,94  ,P5,  P'.  ,  [  F  (  1  A,  I  R  )  ,  t  K  (  I  A,  I  8  I  ,  TC01  Ml  ,  ERR  I  *  I  , 
1,82 
IFISENSE  SWITCH  3)  9S, 10? 
F0RMAT1F8.2,E9.2,F8.4) 
PRINT  9S,P4,P5,B2 

GO  TO  102 

00    9  74    1= I.NM 
00    972    J=1,NHP 
UHCI 1, JI-UHI  1,  J) 
00    9711    L=  I.NCP 
UCCI I.LI-UCI I.LI 


103 

104 

105 
106 


IFISENSE  SWITCH  3)  HI,  132 
F0R*AT</4F8.2,E8.2,F8.2,  3X,  I  3/1  I5F7.  ll )  I 
WRITE  OUTPUT  TAPE  3, 1 3 1 , P 1 , P2 , P3 , P4, P5 , P6, I C , 
IIUHCIJ,  l),UCCINMP-J,NCP) , J  =  1  ,NMI 
P6»P6-P61N 
IB=I8* I 

IFIP6-P6FII  103,103,51 
P6=S6 
IB- I 


P5=P5«P5FA*P51N 

IA=IA.I 

1FIP5FI-P51     105,50,50 

P5=S5 

IFIIC-I)     106,106,108 

IC  =  2 


108 
109 


NMpNM2 

NH=NH2 

NC-NC2 

00     107    K=I,NM 

I 'Mir  ,  I  J     ,ih.  I  »   I 

GO    TO    9 

IFIP4F1-P4)     113, 

P4=P4*P4  IN 


NC=NC-NCI 
UO     III    K=I,NM 
UHIK,  1  l-UHUKI 
GO    TO    8 


1       PAUSE    50 

IFISENSE    SWITCH     I  )     I  I  li ,  I 
I  lit       00     1  15    L=I,NM 
115      UHlLt 1)>UH2(L) 

GO    TO    5 

END 

END 


TO 


UNCLASSIFIED 


UNCLASSIFIED 


